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ABSTRACT 

We present our XMM-Newton observation of the fastest rotating spiral galaxy 
UGC 12591. We detect hot gas halo emission out to 110 kpc from the galaxy 
center, and constrain the halo gas mass to be smaller than 3.5 x 1O 11 M . We also 
measure the temperature of the hot gas as T = 0.64 ± 0.03 keV. Combining our 
X-ray constraints and the near-infrared and radio measurements in the literature, 
we find a baryon mass fraction of 0.03-0.04 in UGC 12591, suggesting a missing 
baryon mass of 75% compared with the cosmo logical mean value. Combined 
with another recent measurement in NGC 1961, the result strongly argues that 
the majority of missing baryons in spiral galaxies does not reside in their hot 
halos. We also find that UGC 12591 lies significantly below the baryonic Tully- 
Fisher relationship. Finally, we find that the baryon fractions of massive spiral 
galaxies are similar to those of galaxy groups with similar masses, indicating that 
the baryon loss is ultimately controlled by the gravitational potential well. The 
cooling radius of this gas halo is small, similar to NGC 1961, which argues that 
the majority of stellar mass of this galaxy is not assembled as a result of cooling 
of this gas halo. 

Subject headings: galaxies: halos — galaxies: individual (UGC 12591) — X-rays: 
galaxies 



Introduction 



Observations show that nearby galaxies are missing most of their baryons (e.g., Hoekstra 
et al. 2005; Heymans et al. 2006; Mandelbaum et al. 2006; Gavazzi et al. 2007; Jiang & 
Kochanek 2007; Bregman 2007) when compared to the cosmological baryon to matter ratio 
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(e.g., f b = 0.171 ± 0.009 from WMAP, Dunkley et al. 2009). For example, the Milky Way 
is missing two thirds of its baryon allotment (e.g., Sakamoto et al. 2003) and less massive 
galaxies have retained less than 10% of their baryons (e.g., Corbelli 2003; Walker et al. 2007). 
This situation is confirmed in other galaxies through a variety of methods (e.g., Hoekstra 
et al. 2005; McGaugh 2007). However, most of these studies do not include the baryon 
mass in the hot gas halo of galaxies, and it is possible that the majority of the missing 
baryons in galaxies actually resides in their hot has halos based on theoretical predictions 
(e.g., White & Frenk 1991; Sommer-Larsen 2006; Fukugita & Peebles 2006; also see the 
review by Benson 2010). This component can be difficult to detect for spiral galaxies due to 
its faintness, especially if the gas density profile is flat. As another possibility, the missing 
baryons from galaxies may have escaped from the potential wells of the galaxies but reside 
in their parent groups or clusters (e.g., Humphrey et al. 2011). Finally, the missing baryons 
can be in the form of warm- hot intergalactic medium (Cen & Okstiker 1999; 2006). Deep 
X-ray observations are needed to distinguish these different scenarios. 

The situation is different for rich galaxy clusters, where the gas mass dominates the 
baryon content, and the baryon fractions in these massive systems are close to the cosmo- 
logical value after combining the gas and stellar baryon contributions (e.g., Vikhlinin et al. 
2006; Allen et al. 2008). The measurements of the baryon fraction in different systems sug- 
gest that the fraction depends on the dynamical mass of the systems: rich clusters retain 
their cosmological allotment of baryons, while galaxies are baryon-poor. We summarize the 
situation in Dai et al. (2010) by combining the archival data points reported in the litera- 
ture and our stacking analysis result using the ROSAT All-Sky Survey data of 4,000 nearby 
galaxy groups and clusters (Dai et al. 2007). We find that the baryon fractions from dwarf 
galaxies to rich galaxy clusters can be fit by a broken power-law model with the break at the 
circular velocity of V c ~ 440 km s -1 . The scatter of the fractions about the mean relation is 
small considering the huge dynamic range of the systems. Further examining the relation, we 
find that the baryon fractions are similar for different systems with similar total masses but 
different compositions. For example, the baryon fractions of poor galaxy groups, where the 
baryon mass is still dominated by the gas mass, are close to those of massive galaxies, where 
the baryon mass is dominated by the stellar mass. Such a coincidence is puzzling considering 
the differences between their mass compositions and energy feedback mechanisms. 

To test whether the missing baryons reside in galaxy halos and further constrain the 
baryon loss in different mass scales, we focus on the massive galaxy UGC 12591. UGC 12591 
is a spiral galaxy with the largest measured rotational velocity to date (466-500 km/sec, 
Giovanelli et al. 1986; Paturel 2003). To appreciate this galaxy, its optical-IR luminosity is 
nine times that of M31. Nine big spirals were in a single group is richer than Hickson groups 
and is about half of the Fornax cluster, but much more compact. In this paper, we combine 
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our XMM-Newton observation of UGC 12591 with 2MASS and other data to determine the 
baryon fraction and the composition of this massive spiral galaxy. Throughout the paper, 
we use the cosmological parameters from WMAP with H = 72 km s _1 Mpc -1 , Q m = 0.26, 
and Vl A = 0.74 (Dunkley et al. 2009). 



2. XMM-Newton Observation and Data Reduction 

We observed UGC 12591 with the XMM-Newton X-ray Observatory (Jansen et al. 2001) 
on December 15, 2008 for 80 ks. The observation log is listed in Table [TJ We reprocessed the 
PN and MOS data using the SAS9 .0.0 software tools epchain and emchain, and filtered the 
events with the patterns < 4 and < 12 for the PN and MOS chips, respectively. We filtered 
background flares by excluding the intervals with background count rates of CR > 0.4 cnt s _1 
and CR > 0.2 cnt s -1 in the 10-12 keV band in PN and MOS observations, respectively, 
following the standard suggestion from SAS. We also applied a low energy flare filter to 
exclude flares in the 0.6-1.4 keV band (Tabled]). We obtained net exposure times of 31.2, 
50.0 and 46.0 ks for PN, MOS1 and MOS2 CCDs, respectively. We detected serendipitous 
sources in the field using the SAS tools, and masked the serendipitous source regions in the 
subsequent analysis. The central region of UGC 12591 is clearly detected in X-rays in all 
PN and MOS CCDs. However, a detailed analysis of the background is needed to determine 
the extent to which the outer region is detected. 



2.1. Background Determination 

We estimated the PN and MOS backgrounds by directly fitting the background spec- 
tra (e.g., Kuntz & Snowden 2008; Leccardi & Molendi 2008). Here, we briefly describe 
the method. The XMM-Newton background is composed of several components induced by 
soft protons, cosmic rays, galactic, and extra-galactic background photons (e.g., Carter & 



Table 1. XMM-Newton Observation of UGC 12591 



Sequence 


Observation 


CCD 


Total 


Effective 


Soft Flare 


Hard Flare 


Number 


Time 




Exposure 


Exposure 


Filters 


Filters 








(ks) 


(ks) 


(cnt s" 1 ) 


(cnt s _1 ) 


0553870101 


2008-12-15 


PN 


80 


31.2 


> 0.9 in 0.6-1.4 keV 


> 0.4 in 10-12 keV 






MOS1 


80 


50.0 


> 0.5 in 0.6-1.4 keV 


> 0.2 in 10-12 keV 






MOS2 


80 


46.0 


> 0.5 in 0.6-1.4 keV 


> 0.2 in 10-12 keV 
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Read 2007). After filtering the background flares in the hard and the soft energy bands, we 
had removed most of the variable soft proton background. We then extracted the spectrum 
from an outer annulus background region, and the remaining background components in 
this spectrum is from quiescent soft protons, cosmic rays, galactic, and extra-galactic back- 
ground photons. The cosmic ray background can be subtracted by comparing the surface 
brightnesses in the outer ring and out-of-FOV (OFOV) regions. We can then decompose the 
contributions from the photon background and quiescent soft proton background through 
spectral analysis. We fixed the spectral index of the extra-galactic background as T = 1.4 
and Galactic background with T = 0.197 keV, and allowed the normalizations of these two 
components to vary during the spectral fits. Figure [1] shows the PN and MOS spectra of the 
outer ring background region, where we used the OFOV spectra as backgrounds. The spectra 
are well fit by the combination of photon background, quiescent soft proton background, and 
residual metal emission lines. We estimated that the ratios between the quiescent soft proton 
background and photon background are 0.676, 0.200, and 0.087 for PN, MOS1, and MOS2, 
respectively, in the 0.6-1.4 keV band. After this step, we can scale the two components to 
inner regions using the corresponding vignetting profiles. 



2.2. Surface Brightness Profile of UGC 12591 

We measured the surface brightness profile of UGC 12591 by extracting events in con- 
secutive annuli of 0.'5 widths (13.6 kpc). For the region within 1' from the center, we further 
divided the annuli using 0.'25 widths. We measured average exposure times within each an- 
nulus. To increase the signal-to-noise ratio of the surface brightness profile, we measured 
it in the soft X-ray (0.6-1.4 keV) band. We first subtracted the surface brightness of the 




Fig. 1. — XMM-Newton PN and MOS spectra extracted from an outer annulus background 
region. We fit the spectra using the OFOV spectra as backgrounds, and modeled the spectra 
as a combination of photon, quiescent soft proton, and residual metal emission line induced 
backgrounds. 
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Fig. 2.— Surface brightness profile of UGC 12591 from the XMM-Newton PN observation. 
The circles represent the raw surface brightness profile, and the triangle represents the OFOV 
surface brightness. The cross symbols are results of subtracting the raw surface brightness 
profile by the OFOV surface brightness. The dashed line is our model for the background 
including the contributions of photon, quiescent soft proton, and residual metal emission 
lines, which fit the outer regions well. The squares and their associated error-bars denote 
the net surface brightness profile and its uncertainty of UGC 12591 (crosses minus the dashed 
line) corrected for vignetting. We also plot the uncertainties as dotted lines. 
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Fig. 3. — Surface brightness profile of UGC 12591 from the XMM-Newton MOS1 observation. 
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Fig. 4. — Surface brightness profile of UGC 12591 from the XMM-Newton MOS2 observation. 



-8 - 



OFOV region, which is un-vignetted, from the surface brightness profile. We then subtracted 
the photon and non-photon induced background based on our background decompositions 
obtained by fitting the outer background spectra. We plot the net surface brightness profiles 
and their uncertainties of UGC 12591 after correcting for vignetting measured from PN and 
MOS CCDs, in Figures HH respectively. UGC 12591 is detected out to ~3' in the PN data 
and ~1.'5 in the MOS data. We also extracted the surface brightness profile of UGC 12591 
in the hard X-ray band between 2-8 keV to measure the contribution from point sources. 
The hard component from point sources is detected to 1' and 0.'5, respectively, in the PN 
and MOS data. After scaling the 2-8 keV count rate to the 0.6-1.4 keV count rate for the 
point source contribution, assuming T = 1.56 (Irwin et al. 2003), we found the contribution 
from point sources in the 0.6-1.4 keV band flux is 45% within the central V region. 



2.3. Spectral Analysis 

We extracted the central spectra of UGC 12591 within 50" from the PN and MOS data. 
The background regions were chosen in a region with large off-axis angles to avoid possible 
contamination from the extended halo emission detected in UGC 12591. We built the rmf 
and arf files separately for the source and background spectra. We modeled the source 
spectra by combining an APEC model for the extended gas emission and a power-law model 
for point sources, modified by Galactic absorption and absorption at the source. We fixed 
the photon index of the power-law component as V = 1.56 and scaled its normalization 
based on the 2-8 keV count rate of the central 50" region, since the expected gas emission 
(the APEC model) in thi s band is negligible. We als o fixed the Galactic absorption based 



on the value obtained by iDickey fc Lockman I (119901 ) . Since our spectra cannot constrain 



the gas metallicity reliably, we fixed it at 0.5 solar metallicity. The free parameters left are 
the temperature and normalization of the gas emission and the Ah column density of the 
absorption at the source. Figure [5] shows our simultaneous fit to the PN and MOS spectra 
using this model. We obtained accepted fit to the spectra with Ax 2 /dof = 131.5/115, and 
constrained the gas temperature as T = 0.64 ± 0.03 keV and the absorption at the source 
as Ah = (5 ± 2) x 10 20 cm -2 . The 0.2-10 keV luminosity for the point source component is 
1.1 x 10 41 erg s _1 . At this luminosity, the central source is possibly dominated by an AGN. 



3. Mass of the Hot Gaseous Halo 



We fit the surface brightness profile in the 0.6-1.4 keV band using three components, 
the hot gas, AGN/binaries, and the contribution from the stellar population (including 
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Fig. 5.— XMM-Newton PN and MOS spectra of the central region of UGC 12591. We 
fit the spectra using a combination of a power law and an APEC model both modified by 
Galactic absorption and absorption at the source. 
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cataclysmic variables). The AGN/binary contribution can be measured using the surface 
brightness profiles from the 2-8 keV band, as described above. We estimated the X-ray 
emission from the stellar population (excluding the binaries) using the stellar mass to X- 
ray luminosity conversion factor of Revnivtsev et al. (2008). This relation is calibrated for 
old stellar populations, and since UGC 12591 is fairly bulge-dominated, the relation should 
still work approximately. To compute the radial surface brightness profile from the stellar 
population in the X-ray band, we derived a K-band radial surface brightness profile from 
the K-band magnitudes within circles of different angular sizes for this galaxy listed in the 
2MASS Extended Source Catalog. We used a distance modulus of 35.0 and a K-band mass- 
to-light ratio of 0.6 (Bell and de Jong 2001) for this galaxy, and then applied the Revnivtsev 
et al. (2008) conversion to derive an X-ray surface brightness profile. The remaining emission 
attributes to the hot gas in the halo. We fit the surface brightness profile of this emission 
using a standard (3 model. 

Figure [6] shows the XMM-Newton PN and MOS data, after vignetting correction and 
background subtraction, as well as the estimates of the various components we fit to the 
surface brightness profile. It is clear that out to 30 or 40 kpc the X-ray binary and stellar 
components are insufficient to account for all the emission, and therefore that a hot halo 
component is necessary. We only allowed /9-models for the hot gas if they cannot be excluded 
at greater than 95% confidence. This results in a narrow range of acceptable fits for the hot 
halo profile. The true uncertainties in the surface brightness profile are probably somewhat 
larger than the statistical errors, however, due to inevitable systematic errors in the flat 
fielding and background subtraction. This means that the range of formally acceptable fits 
for the hot halo might be a little wider than we have indicated Figure [6j Deviations from 
a simple /3-model for the hot gas distribution are also possible, and we tested by adding a 
second flatter /3 model component in our fits. 

The acceptable fit with the highest enclosed mass has (3 = 0.43, r = 1.0 kpc, and 
Sq = 45.3 counts ksec -1 arcmin -2 . This corresponds to a count rate within a projected radius 
of 50 kpc of 0.0032 counts s _1 . Assuming an APEC model with kT = 0.64, Z = O.5Z , 
and Nu = 5 x 10 20 cm -2 , we constrained the mass within 50 kpc as 4.4 x 1O 9 M and the 
unabsorbed luminosity (0.6-1.4 keV) of 2.3 x 10 40 erg s -1 given a distance of 100 Mpc. The 
acceptable fit with the lowest enclosed mass has (3 = 0.69, r = 8.6 kpc, and 5*0 = 6.7 
counts ksec -1 arcmin -2 . This corresponds to a count rate within a projected radius of 50 
kpc of 0.0028 counts s -1 , yielding a mass of 3.9 x 1O 9 M and an unabsorbed luminosity of 
2.3 x 10 40 erg s -1 . The best-fit profile has (3 = 0.52, r = 4.05 kpc, and So = 11.2 counts 
ksec -1 arcmin -2 . The \ 2 is 4.6 for 6 degrees of freedom. If we integrate these profiles out to 
500 kpc, the fit with the highest enclosed mass contains 2.2 x lO n M with an unabsorbed 
luminosity of 3.8 x 10 40 erg s -1 . The fit with the lowest enclosed mass contains 4.4 x 1O 1O M 
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Fig. 6. — Background-subtracted and un-vignetted XMM-Newton PN and MOS radial sur- 
face brightness profiles of the 0.6-1.4 keV emission around UGC 12591 (the open circles 
denote negative values). The PN emission has been multiplied by a factor of 0.26 to convert 
it into effective counts for the MOS detectors. We fit the surface brightness profile with the 
sum of three components, the hot gas halo, AGN/binaries, and the stellar contribution. The 
black dashed line is the emission from coronally active stars and cataclysmic variables in the 
galaxy, estimated by scaling the 2MASS K-band luminosity. The colored dashed lines are the 
emission from X-ray binaries/AGN, estimated by scaling from the 2-8 keV surface brightness 
profile. Finally, the shaded region denotes the emission from hot gas, which makes up the 
rest of the observed surface brightness profile. The black solid line is the best-fit /3-model 
for the hot halo component with a \ 2 °f 61.1 for 45 degrees of freedom, and the green solid 
line shows our best fit by adding an additional flattered f3 profile. 
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with an unabsorbed luminosity of 3.2 x 10 erg s . 

We also examined the possibility of a higher-entropy halo as predicted by many simu- 
lations (Mailer and Bullock 2004; Kaufmann et al. 2009; Crain et al. 2010; Guedes et al. 
2011). Such a halo would have a flatter density profile than the (3 ~ 0.5 models we find 
above, and it would therefore be both more difficult to detect in emission and also more 
massive than a lower-entropy halo. As in Anderson and Bregman (2011), we chose to model 
a flattened profile with a two-component fit to the data. The "flattened" component is a 
/3-model with fixed (3 = 0.35 and r = 50 kpc, but with free normalization, and the other 
component is a more concentrated /3-model with all three parameters free, used to model 
the emission at smaller radii. For this galaxy, there is very little statistical space left to add 
a flattened profile, so the most mass that can be included by adding a flattened component 
is 3.5 x 1O U M . As before, however, we caution that this statistical constraint depends on 
understanding all the systematic uncertainties perfectly, especially the background in XMM- 
Newton PN/MOS CCDs. These profiles can be independently tested in absorption profiles 
instead of emission due with the linear dependence on density in absorption. In addition, 
we note that the metallicity assumed can be another uncertainty in our analysis since the 
total mass will depend on the metallicity. 



4. Discussion 

4.1. Baryon Mass Components in UGC 12591 

We list the various baryon mass components in UGC 12591 in Table [2j We have con- 
strained the hot gas mass of UGC 12591 using the XMM-Newton observation as (4.1 ±0.3) x 
10 9 M Q within 50 kpc with an average temperature of T = 0.64 ± 0.03 keV. We have also 
constrained the hot gas mass of 1.3 x 1O 11 M within 500 kpc regions using our best fit (3 
model, and the hot gas mass is below 3.5 x lO n M within 500 kpc even if we add another 
flatter (3 model component in our fits. Beside the hot gas mass, there are other baryon mass 
components in the galaxy including the ste ll ar ma ss and cold gas mass components. For the 
cold gas mass component, iGiovanelli et al.l ( 119861 ) measure the H I mass of 5.3 x lO 9 /if 2 2 M 



from the radio data. Assuming the H I mass is 75% of the total cold gas mass, we find that 
the total cold gas mass is M cg = 7.1 x 1O 9 M , where we ignore the contribution from the 
molecular gas component. We estimate the stellar mass of UGC 12591 as (4.5±1.0) x lO n M 
within 29 kpc radius using its K band total magnitude (K = 8.89 mag) from the 2MASS 
Extended Source Catalog and a range of mass-to-light ratio from 0.6 to 0.95 (e.g., Bell et al. 
2003). The 2MASS team calculates the total magnitude by integrating the surface brightness 
profile out to ~ 4 disk scale lengths from the isophotal aperture well below the la noise level. 
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For the mass-to-light ratio, Bell et al. (2003) measure a value of 0.95 as the cosmic mean 
value. However, since UGC 12591 is a late-type galaxy and could have a lower mass-to-light 
ratio, we choose to use 0.78 ± 0.18 in our calculation. We find the largest uncertainties in 
the baryon mass are from the systematical uncertainty in the stellar mass-to-light ratio and 
the flattened gas halo. Combining the two effects, we find an uncertainty of 2.4 x 1O 11 M 
for the total baryon mass within 500 kpc. 

In the central region within ~ 50 kpc, the baryon mass is clearly dominated by the 
stellar mass, and the stellar-to-gas mass ratio is r sg ~ 39. Using the rotational velocity 
of 466-500 km s _1 , we measure a total mass of mtot = 2.7 x 10 12 M Q within 50 kpc and a 
baryon mass fraction of /& ~ 0.17, consistent with the cosmological baryon fraction. Out 
to the 500 kpc region (r 2 oo— 550 kpc), we use the gravitational mass m tot = 1.9 x 1O 13 M 
estimated from the X-ray data, because the rotational curve is only constrained within 



28h 7 2 kpc ( Giovanelli et al.lll986l ). smaller than the total mass m to t = 2. 7 xlO 13 M estimated 
using the rotational curve. The baryon mass within 500 kpc is nib = 5.9 x 10 11 M Q , and we 
measure a baryon fraction of /& ~ 0.03. Considering a second flattened gas component, the 
baryon fraction within 500 kpc can reach to ft, <0.04. Since we use the smaller total mass 
estimate in the calculation, the baryon fraction quoted should be treated as a conservative 
upper limit. The stellar mass component is still more important with r sg = 3.3 within 
500 kpc, or r sg > 1.3 with the additional flatter gas component. 

To summarize, combining our XMM-Newton observation and the 2MASS and radio 
data in the literature, we have constrained that UGC 12591 has lost at least 75% of the 
baryons compared to the cosmological value. The missing baryons do not reside in the hot 
halos for spiral galaxies. Our result confirms the recent measurements in another giant spiral 
NGC 1961 using Chandra by Anderson & Bregman (2011), who find that NGC 1961 has 
also lost 75% of its baryon content. 



4.2. Baryonic Tully-Fisher Relationship 

Using the fastest rotating galaxy UGC 12591, we are able to extend the baryonic Tully- 
Fisher relationship (BTF), a correlation between the baryon mass and rotational velocity of 
galaxies (McGaugh 2005; 2011) to the high rotational velocity regime of 500 km s _1 . We 
plot UGC 12591 in the BTF diagram together with the galaxies in McGaugh (2005) and the 
other massive galaxy NGC1961 (Anderson & Bregman 2011) with a rotational velocity of 
402 km s" 1 in Fi gure [71 Anderson & Bregman (2011) find that NGC1961 deviates slightly 
from the linear fits to BTF. However, the authors caution that the offset can be caused 
by systematic uncertainties. Indeed, assuming a K band mass-to-light of 0.95, NGC 1961 
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Table 2. Gravitational and Baryon Mass Components in UGC 12591 



Component 


within 50 kpc radius 


within 500 kpc radius 


Stellar 


45 ± 10 


45 ±10 


Cold Gas 


0.7 


0.7 


Hot Gas 


0.41 ±0.03 


13 


Including Flattened Hot Gas 


0.55 


< 35 


Gravitational 


270 


1900 


Stellar-to-gas ratio r ag 


39 


3.3 


r sg including the flattened gas 


33 


3.3-1.3 


Baryon Fraction fi 


0.17 


0.03 


/(, including the flattened gas 


0.17 


0.03-0.04 



Note. — The masses are in the unit of 10 Mq. The T200 is at ~ 550 kpc. 
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Fig. 7. — The baryonic Tully-Fisher relationship of McGaugh (2005) adding NGC 1961 
(Anderson & Bregman 2011) and UGC 12591. 
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would be on the BTF relation. Here, UGC 12591 provides another challenge to the linear 
BTF relation, which predicts a baryon mass of 2.2 x 10 12 M o for UGC 12591, whereas we 
measure a baryon mass in the range of 5.9-8.1 xlO n M with an uncertainty of 2.4 x 1O U M . 
Thus, UGC 12591 is 6a below the BTF relation. If the offset from the BTF relation is 
caused by the uncertainties in the K band mass-to-light ratio, a ratio of 3.1 is needed to 
put UGC 12591 on the BTF relation, which is extremely unlikely (e.g., Bell et al. 2003). 
Thus, it is possible that the BTF relation turns over for massive galaxies with v c > 400km s -1 . 
However, measurements from a larger sample of massive spiral galaxies are needed to confirm 
this result. 



4.3. Overall Relationship of Baryon Fractions with Total Mass 

We plot the baryon fractions of UGC 12591 and NGC 1961 against their rotational 
velocities, a proxy for the depth of the gravitational potential well, in Figure El We include 
the data composed by Dai et al. (2010) including the stacked results from optically-selected 
clusters (Dai et al. 2007; 2010), individual galaxy clusters (Vikhlinin et al. 2006; Sun et al. 
2008) corrected for the stellar component (Dai et al. 2010), individual galaxies (McGaugh 
2005; Walker et al. 2007; Stark et al. 2009), elliptical galaxies (Gavazzi et al. 2007) and 
the Milky Way (Sakamoto et al. 2003; Flynn et al. 2006), and the recent addition of gas- 
rich late-type galaxies (Begun et al. 2008; Trachternach et al. 2009) composed by McGaugh 
(2011). This enables us to compare the baryon loss across a large range of systems from 
dwarf galaxies to rich galaxy clusters. For massive spiral galaxies such as UGC 12591 and 
NGC 1961, their dark matter halo masses are close to that of a medium galaxy group, and 
we find that their baryon losses are comparable to the stacked results for galaxy groups in 
Figure El For different systems as spiral galaxies and galaxy groups, the consistency between 
their baryon fractions suggests that the baryon loss is ultimately controlled by the potential 
well of the dark matter halo. However, currently we still lack constraints from individual 
galaxy groups to confirm the stacking results of Dai et al. (2010). 

The overall baryon fractions for all systems can be fit by a broken power law model 
(Dai et al. 2010). With the addition of new measurements, especially from those gas-rich 
late-type galaxies (Begun et al. 2008; Trachternach et al. 2009), we can better fit the power 
law slope for baryon fractions in less massive systems. Thus, we re-fit the data with a broken 
power law model to find that 

0.16(w c /643 km/s) a 
b ~ (l + (v c /QA3 km/s)f /c ' 

where a = 1.26, b = 1.24, and c = 2. The baryon fraction, scales as fb oc -y" _6=002 
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Fig. 8. — The baryon fraction as a function of the gravitational potential well indicated 
by the circular velocity at r2oo- We plot the new measurements from massive spiral galax- 
ies UGC 12591 and NGC 1961 (Anderson & Bregman 2011) over the archival data from 
Sakamoto et al. (2003), McGaugh (2005), Flynn et al. (2006), Vikhlinin et al. (2006), Gavazzi 
et al. (2007), Walker et al. (2007), Stark et al. (2009), Sun et al. (2009), Dai et al. (2010). 
The dotted line is the cosmological baryon fraction measured from CMB, and the dashed 
line is our best fit broken power law model for baryon losses. For massive spiral galaxies like 
UGC 12591 and NGC 1961, the baryon loss is at least 75%. 
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above the break and ex v% =lm26 below the break, and the parameter c in the equation is 
the smoothness of the broken power law model, which is fixed in our fit. Comparing with 
the fit in Dai et al. (2010), we find the major difference lies in the shallower slope fb oc f 1 ' 26 
for the baryon loss in galaxies. We also find a larger break location and a shallower slope 
for galaxy clusters. 



4.4. Cooling of the Gas Halo 

The gas halo is predicted to play an important role in galaxy formation and evolution. 
With our detection of the gas halo emission in UGC 12591, we can estimate the cooling time 
of this hot halo and the implied accretion rate onto the galaxy, which can provide constraints 
on the gas available for new star formation. We define the cooling radius as the radius where 
the cooling time is 10 Gyr, using the expression of the cooling time (Fukugita & Peebles 
2006), 

l.hnkT ^ 1.5kT x 1.92 

T ^ ~~ An e (n - n e ) ~ An e x 0.92 ' ^ ' 

where the latter expression assumes a primeval He abundance resulting in a total particle 
density of n = 1.92n e . For T = 10 6 - 85 K, Z/Z Q = 0.5, and A = 10" 22 - 85 erg cm 3 s" 1 
(Sutherland & Dopita 1993), the cooling radius is at n e = 6.8 x 10~ 4 cm" 3 . For the range 
of best-fit /3-model profiles constrained in this paper, this corresponds to a cooling radius 
between 15.6 and 18.0 kpc, and a hot halo mass of 6.2 — 9.2 x 1O 8 M within that radius. We 
can roughly estimate the cooling time and rate by dividing the thermal energy in the hot gas 
within the cooling radius by the luminosity within that radius, and this yields a wide range 
in cooling time of 2.8-6.3 Gyr for material within the cooling radius, but a fairly narrow 
range in the effective cooling rate of 0.15 — 0.21M Q year" 1 . This halo accretion rate is two 
orders of magnitude too low to assemble the stellar mass of this galaxy within a Hubble time. 
Therefore, significant accretion must have occurred via some other mode, such as cold flows 
or mergers, to produce the stellar mass seen in this galaxy today, confirming the conclusion 
drew in NGC 1961 (Anderson & Bregman 2011). 
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